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jSUMMARY 


The stability of three-dimensional rotating disk flow is 
investigated, including the effects of Coriolis forces and 
streamline curvature. The numerical results show that the 
critical Reynolds number for establishment of stationary vortex 
flow is = 287. These vortices spiral outwards at an angle 
of about 11.2° and transition to turbulence occurs when their 
total amplification is about e^^. It is shown that our anal- 
ysis gives growth rates that compare much better with the avail- 
able experimental results than do results obtained using the 
Orr-Sommerfeld equation. The experimental results tend to 
support the numerical prediction that the number of stationary 
vortices varies directly with the Reynolds number, Our cal- 
culations also indicate the existence of weakly unstable pro- 
pagating (Type II) modes at low Reynolds numbers (R^ - 49) . 


i . INTRODUCTION 

The prediction of transition in three-dimensional boundary 
layers (1-3] is a rsubject of both fundamental and practical im- 
portance in fluid mechanics. Practical interest in the subject 
centers on the design of laminar flow control (LFC) wings that 
promise significant improvement in airplane fuel efficiency. 

At present, the most useful tool for transition prediction in 

N 

such flows is the so-called e method [4], Hefner and Bushnell 
tS] and Malik and Orsssag [6] show that the exponent N (called 
the N factor) is of the order 7-11 when transition occurs 
on LFC swept wings. 

The instability mechanism exhibited in the leading edge 
region of a swept wing is similar to that found in the boundary 
layer on a rotating disk, since both have mean cross-flow pro- 
files with inflection points. More details on the similarities 
between the two flows is given in Ref. [7]. The rotating disk 
flow is amenable to stability analysis in view of von Karman's 
exact steady solution [8] of the Navier-Stokes equations free 
of boundary-layer assumptions. 

Using hot wire techniques, Smith [ 9 ] observed that sinus- 
oidal disturbances appear in a rotating disk boundary layer at 
sufficiently large Reynolds number. About 32 oscillations were 
observed within a disk rotation period and analysis indicated 
that the disturbances propagate at an angle of about 14° relative 
to the outward drawn radius (where the direction of disk ro- 
tation defines positive angles). Later,. in a remarkable study 
using the china-clay technique, Gregory et al [10] observed 


1 


about 30 vortices over the disk spiraling outwards at an angle 
of about 14° (that is, their normals make an angle of about 14° 
with the outward drawn radius). These vortices, which appeared 
stationary relative to the rotating frame of the disk, were first 
observed at a Reynolds number Rg«430 [where R is defined after 
(10) below] < Transition to turbulence occurred at R^«530, The 
stationary vortex flow established due to a rotating disk has 
later been studied by several investigators [11-13], There 
exists, among these studies, considerable controversy over the 
value of the critical Reynolds number which in our view can be 
attributed to the measurement techniques used. There is also 
some confusion over the number of stationary vortices. Fedorov 
et al [13], using visual (Napthalene) and acoustic techniques, 
observed 27-30 vortices at Reynolds numbers R > 387, However, 
at low Reynolds numbers, they observed only 14-16 vortices with 
normals lying at an angle of about 20°. 

Stuart (10] analyzed the linear, inviscid stability of 
rotating disk flow. However, the neglect of viscosity resulted 
in the prediction of 113-140 vortices over the disk which is 
about four times larger than the observed value. Brown [14] 
extended Stuart's work to the viscous case by applying the Orr- 
Sommerfeld equation to disk flow. Using temporal instability 
theory, Brown found R^k 178, which is much less than the observed 
value. Recently, Cebeci and Stewart son [3] solved the Orr- 
Somrrs. 2 rfeld equation for rotating disk profiles using spatial 
stability theory and found R^=176. They, also suggested that 
wave packets propagate in three dimensional flows in such a way 
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that da/d& is real . Using this condition Cebeci and Stewartson 
correlated transition using the e method and found the N 
factor at transition 510) to be about 20 which is much 

higher than that found for LPC wings ([5], 16]). Bushnell 
(private communication) argues that a higher N factor may be 
required for transition in dis); flow than on LPC wings because 
the boundary layer is rotating with the dislc while the external 
disturbances in the surroundings are not. Consequently, there 
is no direct coupling between the external disturbances and the 
instability waves in the rotating disk boundary layer. 

The Orr-Sommerfeld stability equation neglects the effects 
of Coriolis forces, streamwise curvature, and nonparallel flow. 

In Ekman layer flow, Lilly [15] has shown that the Coriolis 
force has a significant effect on stability at low Reynolds 
numbers. Lilly showed that the critical Reynolds number for 
appearance of stationary vortices is higher (R^«115) when the 
Coriolis force is included in the analysis than when it is ne- 
glected (R =*85) . In addition to the stationary vortices, he 
showed the existence of another (Type II) instability caused 
by the Coriolis force at much lower Reynolds numbers. Such an 
instability mechanism was also observed in the Ekman layer 
experiments of Paller and Kay lor [16] and Tatro and Mollo- 
Christensen [17]. The Ekman layer and the rotating disk are 
similar in that both are three-dimensional layer flows in which 
rotation plays a significant role. Lilly's results suggest that 
the inclusion of the Coriolis force in the stability analysis 
of rotating disk flow may also lead to a higher critical Reynolds 
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nun±Mtr for stationary vortices which is in better agreement with 
observations. 

In this paper we present a stability analysis of rotating 
disk flow in which the effects of Coriolis force and streamline 
curvature are included. The resulting sixth-order system is 
solved numerically by a Chebyshev spectral method (18,191. We 
also follow the evolution of the disturbance modes using the 
envelope method [1,6] and calculate the N factor at transition. 
The work of Kobayashi et ai (12], which appeared during the 
final stages of the present study, also includes the effects 
of the Coriolis force and streamline curvature. We will comment 
on this work in Sec. 6. 

2. THEORETICAL AimLYSIS 

Consider an infinite plane rotating about its axis with 
angular velocity H . We take cylindrical coordinates r*,6,s* 
with 2* ?! 0 being the plane of the disk and assume the fluid 
to lie in the half-space 2*>0. Let p, u, v, w, denote the 
steady state values of pressure and velocity in the r*,0,5!* 
directions, respectively, in the rotating coordinate frame. 

Von Karman's exact solution [6] of the Navier-stokes equations 
for steady laminar rotating disk flow is obtained by setting 

u - r*nr( 2 ), V = r*nG(z), w =/vQH (z) , p = pvS]p(z), (1) 

where z = z * . The Navier-Stokes equations reduce to the. 
following equations for F,G,H and P: 

- (G+1)^ + F'H - F" = 0 ' (2) 

(3) 


2F(G+1) + G' H - G" 


0 
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p' + HH’ - H" * 0 ( 4 ) 

2F + H' » 0 (5) 

where the prime denotes differentiation with respect to z. The 
boundary conditions are 

F » 0, G « 0, H * 0 (z « 0) (6) 

F * 0, G » - 1 (z •+ ») 

Now we study the evolution of infinitesimal small dis- 
turbances imposed on the st^^ady flow governed by Eqs. (1) - (5). 

Let r* be the radial location near which the analysis is to be 
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made. Using r* H as the reference velocity, 6 * “ 


as the 


reference length, and pr| as the reference pressure, the 
instantaneous nondimensional velocities u,v,w and pressure p 
can be written as 

u(r,6,z,t) = I F(z) + u(r,e,z,t) (7) 

v(r,0,z,t) = I G(z) +v(r,0,x,t) (8) 

w(r,0,z,t) = I H(z) + w{r,0,z,t) (9) 


p(r,0 ,z ,t) 


P(z) + p(r,0,z,t) 


Here the nondimensional radius is r = r*/?27v, the Reynolds number 
is R = r* /n/v , and r^ corresponds to r = R. 

0 lU 

Substituting Eqs. (7) - (10) in the Navier-Stokes equations 
and linearizing with respect to the perturbations gives: 

i " I " i " I i ^ I ^ I ^ - I V t I F* ; (11) 


r ^ R ^ 


. 3^u . 1 3u 2 3v u , 

rr r I? " T w ‘ 

3z r r 
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i 




( 12 ) 


1 iE o. 1 . 1 3^v . 9~v . 1 3v , 2 9u V , 

r w •*■ R 717 ■*■ Z7 r 3r w •• T-J 


3^v 


1 3v . 2 3u 


3r' 


36' 


3z‘ 


^ + EFiH + G 3 !l+H^ + l!lw 

8t R 3i‘ R 36 R 32 R 


3z ^ 


,3 w . 1 3^w 

[ — "■ » + —K -IK 

3r r 30 


2~ 

a. 5 W . 

T y T 

32 


i 1 Hi 
r 3r^ 


(X3) 


3r r 30 3z r 


(14) 


The houndary conditions are that u,v, and w vanish at z * (J,». 
For R?>1, the system (11) - (14) may be consistently 

approximated by replacing factors of r by R and neglecting terms 

-2 

of order R and smaller* The replacement of r by R at this 


stage of the calculation implies that we neglect some nonparallel 

-2 

flow effects. The neglect of terms of order R and smaller 
has little effect on the results discussed below, as we verified 
by computations in which they were included. 

Replacing factors of r by R in (11) (14) gives a set of 

equations that are separable in r,0,t so that the perturbation 
quantities may be assumed to have the form 

(u,v,w,p) = (f (z) ,h(z) ,(J)(z) ,TT(z))e^^°‘^+®^®"'^'‘^^ (15) 

With this assumption, Eqs. (11) - (14) become (not yet dropping 

-2 

terms of order R ) 

i(aF+3G-a))f+F'(;) + iaiT = | [f'-X^f - Ff + 2(G+l)h - Hf’) 


+ ^ [iaf - 2i 6 h)- f 
R^ r"* 


(16) 
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( 17 ) 


i(oF + $G-w)h + ifJir » | lh*‘-X^h-Ph-2 (G+1) f-Hh' J 

+ — ^ [Xoih + 2i 6 •* h 

« R 

i(ap + 0G-u) 4» + tt' * “ ^ ♦ (18) 

R 

(ia+|)f ^ igh + 4' *0 (19) 

where +■ &^* 

Eliminating tr from (16) •- (18) by means of (19) gives, 

-2 

neglecting terms of order R and smaller, 

I li(D^-X^) (D^-X? + R(aP+eG-w) (D^-X^) •• R(aF"+0G'') -iHD(D^-X^) 

- iH' (D^-P) - iFD^]4+|CG+l)P + 2GM n »= 0 (20) 

I [2(G+l)D-iR(aG’-0F') ] 4+|li(D^-l^) + R(aF+BG-uj) - iHD - iF]n » 0 

( 21 ) 

and where D = d/dz, a^a-i/R,X^ s= aa+0^ and n-ah - 0f is pro- 
portional to the z -component of the perturbation vorticity. The 
final result (20) - (21) is a consistent set of stability equations 
valid to order R 

The boundary conditions for the sixth order system (20) - 
(21) are 

4(0) » 4'(0) = n(0) 0 , 4 ( 0 ?) = 4M“) = nC") = 0 (22) 

Note that if the Coriolis force and streamline curvature 
effects are neglected, the above system reduces to the fourth- 
order Orr- Sommer f eld equation; 

[i(D^-X^)^ + R(aFteG-w) (D^-X^) - R(aP"+BG")]4 *= 0 (23) 
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In Sftc. 6 1 wc report munericaX rcsulte for both the sixth 
order system (20) - (21) and the fourth order equation (23) in 
order to study the effect of Coriolis force and streamline 
curvature terms on the stability of flow due to a rotating disk. 

3. NUH'iSRlCAL METHOD 

We solve the Orr- Sommer £e Id Equation (23) using the computer 
code SALLY (11 that employs a spectral method based on Chebyshev 
polynomials [18-19] . Here we extend the method to solve the 
sixth order system (20) - (22). 

The boundary layer coordinate z, 0 < z < « is mapped into 
the finite interval -Is 5<1 by the algebraic mapping 


£ B 2 — 1 

^ z+L 

where L is a scale parameter chosen 
of points in C* Then $ (z) and n(z) 
finite Chebyshev polynomial series 


(24) 

to optimize the distribution 
are approximated as the 


M 

$(z) “ S a T (U (25) 

n»0 


M 

n(z) S T„ (C) (26) 

n=0 ^ " 


Substituting (25 ) , (26) # in (20) - (22) and collocating 
[19] at the discrete points Cj = cos (0 1 3 £ M) gives the 

algebraic eigenvalue problem, 




( 27 ) 


\ 
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where A and B are 2(M4-1) x 2(M4>1) matrices* The eigenvalue problem 
(27) is solved globally (if a guess for an eigenvalue is not 
available) by a generalized QF algorithm or locally (if a good 
guess is avai table) by inverse Rayleigh iteration (20). The 
resulting scheme is very efficient and accurate* In the present 
calculations the optimum value of the scaling parameter L was 
found to be about 1.8. In most of the calculations reported 
below, M » 34 so 35 Chfbyshev polynomials were used. 

The accuracy of the method was tested in several ways. 

First, the number of retained polynomials,. M -f 1 was varied to 
check the accuracy of the eigenvalues and eigenfunctions. Then, 
calculations were made for the stability of Bkman flow. Com- 
parisons were made with the results obtained by Lilly [15] 
for R “ 65, 110, 150, 300 and 500 with good agreement. 

For rotating disk flow, the global method gives only one 
unstable eigenvalue (lm(w)>0) for R t 150 that is insensitive 
to M. However, spurious unstable modes appear for lower R 
which are discarded as unphysical because they are very in- 
sensitive to M. 

4. TRANSITION PREDICTION USING THE e^ METHOD 

In three-dimensional flow, the dispersion relation is given 
by the complex relation 

0 ) = w(a,6) (28) 

where a, 8, and w are, in general, complex. Therefore, there are 
four arbitrary real parameters among a, 8, and u. There are several 
ways (6] to remove this arbitrariness. In the present study we 
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employ the envelope method [Hi Here the four conditions ere 
obtained by using temporal stability theory [in which Im(a) » 
Im'tS) * 0] and by maximizing Im(w) with respect to at fixed 
Re (w) . The N factor is then given by 


N 


s„ 

r 


Im(d)) 

r"vv~i 


ds 


(29) 


where is the (complex) group velocity and s is the 

arc length along the curve whose tangent is the real part of the 
group velocity. Noting that 


ds s* 



/(Re(w^))^ + (Re(Wg))^ 


dR 

Re(w^) 


Re ( Vg) 


Eg. (29) can be written as 
Rt 

/ Im(oj) 

Rj, Re(w^) 


N 


(30) 


Here the subscripts C and T indicate critical (linearly 
unstable) and onset of transition# respectively. 


5. RESUETS AND DISCUSSION 

Critical and Transition Reynolds Numbers 

Some of the available experimental data for critical and 
transition Reynolds numbers are given in Table 1. It is apparent 
that there is considerable variation of the observed critical 
and transition points. We believe that the variation can be 
attributed to either the different measurement techniques used 
in the experiments or to the influence of external disturbances. 
Using the Orr-Sommerfeld equation, we obtained a critical Reynolds 


number which is in good agreement with the theoretical 

results of Brown (14] and Cebeci and Stewartson (3] but is con- 
siderably less than the observed values. The value of the critical 
Reynolds number for stationary disturbances is significantly im- 
proved when the effects of Coriolis forces and streamline curvature 
are included. Our calculated critical Reynolds number of 287 
for disturbances of this kind is in excellent agreement with the 
results obtained in [21] and [12] using hot wire techniques. 
Kobayashi et al [12] also pex*formed a theoretical analysis in 
which some of the effects of Coriolis forces and streamline 
curvature were considered. They calculated a critical Reynolds 
number of 261. 

In order to correlate transition using stability theory, 
one has to know the experimental location of the onset of tran- 
sition. The transition Reynolds numbers usually given for ex- 
periments (see Table 1) are the locations where transition is 
complete. Gregory and Walker [22] showed that, for a slotted 
rotating disk, the transition region is composed of two sub- 
regions: (i) a vortex region and (ii) an intermediate turbulent 

region where the intermittancy factor y varies from 0 to 1. 
Stability theory is only applicable up to the point where the 
first turbulence burst appears (y = 0). Gregory and Walker 
obtained R = 505 and 524 for y = 0 and y = 1, respectively. 

Chin and Litt [23], using an electrochemical technique, observed 
that the transition was complete at R = 592. They also observed 
that the vortices start breaking down into turbulence at = 510. 
We believe that this result should be taken as the relevant 
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location for the onset of transition for the purposes of comparison 
with stability theory, further evidence that the initiation of 
transition occurs at « 510 is provided by Kobayashi et al 
[12] who observe that the disturbances are non-linear at R « 500. 
Usually in boundary layer flows the non-linear region is relatively 
narrow so the onset of transition soon follows. Further, 

Pederov et al [13] observed turbulent flow at R *= 515. In the 
experiments of Ref. [21] turbulent spots first appeared at R = 
513-526. On the basis of the evidence provided by these in- 
vestigators [12,13,21-23] we take R^ = 513 as the location of 
the onset of transition for the purpose of applying the e^ 
correlation method. 

Growth of Infinitesimal Disturbances 

Disturbances of all frequencies may be present in natural 
transition. We follow the evolution of several different mo(” 
and the one which gives the highest integrated growth factor is 
used to correlate transition. Stationary disturbances were found 
to give the highest N factor for rotating disk flow over all 
positive real frequencies. Disturbances with negative phase 
velocities can give slightly higher N factors but they are of 
no consequence in the process of transition. 

It was shown in [6] that the envelope method is a reliable 
tool for transition prediction in three dimensional flows. 

First, we report calculations using the Orr-Sommerfeld equation. 

The resulting N factors are compared with those of Cebeci and 
Stewartson [3] in Figure 1. It is evident that their method 
predicts N==20 at transition (R = 513) While the present (envelope) 
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method gives N « 22 at transition. Cebeci and Stewartson [3J 
used spatial stability theory and in order to remove arbitrar- 
iness among the parameters of equation (28) , they imposed the 
condition 




(31) 


where » Im(a) , = Re (3). In order to simplify their com- 

putations, Cebeci and Stewartson assumed that the maximum growth 
rate at any R>Rq is independent of the growth direction. We 
believe that had their growth rates been maximized over all 
possible growth directions their N factor at transition would 
be in better agreement with the present predictions using the 
Orr-Sommerfeld equation. 

In Fig. 2 we plot calculated temporal growth rates (Im w) 
for stationary vortices. It can be seen that the inclusion of 
streamline curvature and Coriolis forces have a significant 
stabilizing effect. Calculations with only Coriolis terms (as 
done by Lilly [15] for Ekman flow) were also made. These results 
indicate that streamline curvature effects must also be included 
in order to model properly the physical problem. 

Since the instability is spatial in nature, we transform 
temporal growth rates to spatial growth rates a using the group 
velocity transformation 


T - lm(m) 
Re(u^) 


(32) 


The spatial growth rates are plotted in Figure 3. The 
effect of streamline curvature and Coriolis forces is found to 
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be very aignif leant. 

Integrated growth rate (N factor) results are presented in 
Figure 4. The present stability theory gives N « 10.7 which is 
close to the value N * 9 for two dimensional flows and is in the 
range of values found for swept wings [6). It is apparent that 
there is a very significant effect on the predicted transition 
N factor when the effects of Coriolis forces and streamline 
curvature are included. The resulting N factors are much lower 
than those obtained by conventional stability theory where only 
the Orr- Sommer £e Id equation is solved. 

Also presented in Figure 7 are experimental results from 
Ref. [21] for the factor. The experimental amplification rate 
of the stationary vortices was determined by integrating the 
rms fluctuation spectrum in a narrow band centered on a frequency 
determined by the observed number of vortices i as described in 
[21]. Due to the problem with disk vibrations, the disturbance 
amplitude (Aq) at the minimum critical Reynolds number could 
not be measured with any certainty. It was assumed that A^ = 1 
and the resultant data were shifted at constant Reynolds number 
to match the theoretical growth curve . The data are seen to be 
in a fair agreement with the present theory over the range 400 
<R<500. The significant deviation of the data for R<400 is 
attributed to disk vibrations. The falling off of the data for 
R>500 is due to the highly nonlinear nature of the flow in this 
region and breakdown to turbulence. 

Although shifting the level of the data because of un- 
certainty in the value of Aq can be questioned and the resultant 
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data can not be used to show conclusive support for the theory, 
the fact that the slope of the experimental curve (o * dN/dR) 
matches the present theory in the range 400<R<500 is very 
encouraging i 

Orientation and Number of Vortices 

In the envelope method we maximize growth rates Im(w) for 
stationary vortices over all possible wave angles. We find that 
the vortex spirals make an angle of about 11.2° with the negative 
of the direction of disk rotation. This is within t*ie range 
of the experimental value of 11-14° [10-13]. The predicted angle 
varies very slowly with R. 

It can be shown that the number of vortices is given by 

n = BR (33) 

where B is defined in (15). Our calculations show that B remains 
almost constant (siO.0698) so that n varies nearly linearly with 
R. The numerically predicted variation of number of stationary 
vortices is plotted in Figure 5 along with the experimental data 
[21]. The experimental results show that there are about 21 
vortices at a Reynolds number of 294. The number of vortices in- 
creases radially and there appear about 29-31 vortices at R =: 513. 
This radial increase in n is apparently due to the branching of 
individual vortices with increasing R. Although the data fall 
well below the theoretical curve at high Reynolds numbers, it does 
substantiate the overall trend predicted by the theory. 

Previous investigators [10,12,13,25] have not suggested 
any direct correlation between the number of vortices and the 
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|ieynp:|.4a pi^nbers. A pareful analysis of thplr rlata (Figure 5 ) , 
however, shows that their results do not necessarily disagree 
with the present findings. 

Parallel or Type II Instability 

Lilly [15] presented numerical solutions of the Ekman layer 
problem and included the effect of Coriolis forces in his analysis. 
He found that at very low Reynolds numbers an instability mechanism 
exists whose disturbances are different from the stationary 
disturbances described previously. Lilly called this "parallel 
instability” and suggested that it is of viscous type since it 
vanishes at high Reynolds numbers. He found that the critical 
Reynolds number for these fast moving disturbances is 55 and 
the resulting modes are oriented at small negative angles. The 
orientation angle at the critical point is -23° which decreases 
in magnitude as the Reynolds number increases. A similar in- 
stability mechanism was detected in the experiments of Faller and 
Kaylor [16] (who called it a type II instability) and Tatro 
and Mollo-Christe.nsen [17]. 

In our calculations, we also find travelling disturbances 
(f ~ 100 Hz relative to the disk) at Reynolds numbers much lower 
than the critical Reynolds number for stationary disturbances. 

The critical Reynolds number for travelling disturbances is 
calculated to be about 49. These disturbances appear to have 
characteristics similar to those reported by Lilly [15]. 

6. CONCLUSIONS 

The growth of instabilities in the three-dimensional flow 
due to a rotating disk is studied using linear stability theory 

c)RIGINAL PAGE IS 
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including the effects of Coriolis forces and streamline curvature. 
Using the e^ analysis for transition prediction in these three 
dimensional boundary layers the N value is the order 11. Recent 
experiments [21] tend to support the theoretical prediction that 
the number of stationary vortices increases with R. 
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TABLE 1. EXPERIMENTAL CRITICAL AND TRANSITION REYNOLDS NUMBERS FOR 

STATIONARY VORTEX FLOW 


Reynolds Number 

Onset of Method of 


Investigators 

Critical 

Transition 

Transition 

Investigation 

smith [9] (1947) 

460 

557 

- 

hot-wire probe 

Gregory et al 
(10) (1955) 

430 

530 

- 

visual (China 
clay technique) 

Cobb (t Saunders 
(241 (1956) 

447 

490 

- 

heat transfer 
from the disk 

Gregory & WaliKer 
[22] (1960) 

367 

524 

505 

acoustical 
slotted disk 

Chin & Litt 
(233 (1972) 

412 

592 

510 

mass transfer 
coefficient 
using electro- 
chemical 

Fedorov et al 
(133 (1976) 

387 

515 


visual 

(Napthalene) , 
acoustical 

Clarlcson et al 
[11] (1980) 

532-621 

562-680 

- 

visual 

(dye in water) 

Kobayashi et al 
[12] (1980) 

297 

566 

500 (non- 
linear) 

hot wire probe 

Malik et al (21} 

294 

- 

513-526 

hot wire probe 
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FIGURE a 


FIGURE 3 


FIGURE 4 


Integrated growth factor using Orr-SonunerfeXd 
equation 

(a) Present calculations 

(b) Cebeci and Stewartson f3J 

Temporal growth rates for stationary vortices 

(a) Orr-Sommerfeld equation 

(b) Orr-Sommerfeld equation with Coriolis force 
effects included 

(c| Orr-Sommerfeld equation with Coriolis force 
and streamline curvature effects included 

Spatial growth rates for stationary vortices 

(a) Orr-Sommerfeld equation 

(b) Present theory 

Integrated growth factor for staiionary vortices 

(a) Orr-Sommerfeld equation 

(b) Present theory 

Number of stationary vortices as a function of 

R. 
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